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We present the results of a search for supersymmetry with gauge-mediated breaking and x5 ~* 7^ in the 
77-l-missing transverse energy final state. In 2.6±0.2 fb~^ of pp collisions at y^=1-96 TeV recorded by the 
CDF II detector we observe no candidate events, consistent with a standard model background expectation of 
1.4it0.4 events. We set limits on the cross section at the 95% C.L. and place the world's best limit of 149 GeV/c^ 
on the Y? mass at t-o<SC1 ns. We also exclude regions in the y? mass-lifetime plane for r-o!5~2 ns. 



1. Introduction 

For theoretical reasons and because of the 
'ee77-|-missing transverse energy {frY candidate 
event recorded by the CDF detector in RUN I 2], 
there is a compelhng rationale to search in high en- 
ergy collisions for the production of heavy new parti- 
cles that decay producing the signature of 77 -f 

An example of a theory that would produce such 
events is gauge mediated supersymmetry breaking 
(GMSB) [ij with xi -> iG where the Xi is the light- 
est neutralino and the next-to-lightest supersymmet- 
ric particle (NLSP) and the G is a gravitino which 
is the lightest supersymmetric particle (LSP), giving 
rise to fr by leaving the detector without depositing 
any energy. The G also provides a warm dark matter 
candidate that is both consistent with inflation and 
astronomical observations 0]. 

In these models, above the current limits from re- 
cent experiments [J|, the is restricted to be well 
above 100 GeV and is favored to have a lifetime on 
the order of a nanosecond; the G is restricted to have 
a mass in the range 0.5<mg,<1.5 keV/c^ [EJ. At the 
Tevatron sparticle production is predicted to be pri- 
marily into gaugino pairs, and the Xi mass (m-^o) 
and lifetime (r^o ) are the two most important param- 
eters in determining the final states and their kine- 
matics 1]. Depending on how many of the two Xi's 
decay inside the detector, the event has the signature 
77 -|- fjT, 7 + pr or fjT with one or more additional 
high Et particles from the other gaugino pairs. Dif- 
ferent search strategies are required fory^ lifetimes 
above and below about a nanosecond Q. Previous 
searches have been performed for low lifetime models 
in 77 -h It [i, 01 and nanosecond lifetime models in 
the delayed 7 -I- jet + |t 0, [1] final state. 

In this analysis we focus on the 77 -|- fjr final state, 
as recommended in for low lifetime, high-mass 
models of the Xi- The new features of our analysis 
since the last 77 search with 202 pb~^ using 

the CDF detector jl0| are to use the EMTiming sys- 
tem [ll| and a new Met Resolution Model [l^] ■ We also 
use 13 times the data (2.6 fb~^). These additions sig- 
nificantly enhance our rejection of backgrounds from 



instrumental and non-collision sources, which allows 
us to considerably extend the sensitivity of the search 
for large ^5 masses compared to previous Tevatron 
searches |^ . We also extend the search by considering 
Xi lifetimes up to 2 ns which are favored for larger 

Xi . . 

Our analysis begins by defining a preselection sam- 
ple by selecting events with two isolated, central 
(|?7| < 1.0) photons with Et > 13 GeV. AU candi- 
dates are required to pass the standard CDF dipho- 
ton triggers, global event selection, standard photon 
ID, and non-collision background rejection require- 
ments [i,[i3- 

The final signal region for this analysis is defined 
by the subsample of preselection events that also pass 
a set of optimized final kinematic requirements. The 
methods for determining the background in the sig- 
nal region are based on a combination of data and 
Monte Carlo (MC) simulation and allow for a large va- 
riety of potential final sets of kinematic requirements. 
We perform an a priori analysis in the sense that we 
blind the signal region and select the final event re- 
quirements based on the signal and background expec- 
tations alone. We optimize our predicted sensitivity 
using a simulation of our GMSB model. We then cal- 
culate, for each GMSB parameter point, the lowest, 
expected 95% C.L. cross section limit in the no-signal 
scenario [3] as a function of the following event vari- 
ables: MetSig, A0(7i,72), and Ht, each of which will 
be described in Section [5l 

2. Data Selection 

The analysis is based on 2.59±0.16 fb^^ of data de- 
livered to the CDF detector in Run II. The analysis 
selection begins with events that pass the CDF dipho- 
ton triggers which is effectively 100% efficient for the 
final diphoton selection requirements [l3| ■ We require 
both highest- i?T photons to be in the fiducial part of 
the detector with \ri\ < 1.1, pass the standard photon 
ID and isolation requirements and have E^ > 13 GeV. 
In addition to the standard photon ID requirements 
we have added additional requirements to suppress 
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photomultiplier tube (PMT) high- voltage breakdowns 
("spikes") and electron rejection requirements [Sj] 
to remove events where an electron fakes a prompt 
photon (Phoenix tracking rejection). Each event is 
required to have at least one high quality vertex with 
\zvt\ <60 cm. The Et of all calorimeter objects 
(individual towers, photons, electrons, and jets) are 
calculated with respect to the highest ^ Pt vertex. 
However, an incorrect vertex can be selected when 
two or more collisions occur in one beam-bunch cross- 
ing, making it possible that the highest reconstructed 
^ Pt vertex does not produce the photons. If assign- 
ing the photons to a different vertex lowers the frr, 
we take that frp and the photon Et^s to be from that 
vertex for all calculations (Vertex Re-assignment). 

Additional standard selection requirements are 
placed to reduce non-collision backgrounds, such as 
cosmic rays and beam- related (beam halo) effects Q. 
We also apply frp quality requirements (cleanup) to 
remove events if there is evidence that the second pho- 
ton (72) or a jet is partially lost in a crack between 
detector components Our pre-selection sample 

consists of 38,053 events left after all the quality, ID 
and cleanup requirements are applied [l^ . 



3. Backgrounds 

There are three major sources of background for 
j^+fjT events: QCD events with fake electroweak 
events with real fJr, and non-collision events (PMT 
spikes, cosmic ray or beam-halo events where one or 
more of the photons and fJr are not related to the 
collision) . 

Standard Model QCD sources, 77, 7 — jet 
11 fake, and jet- jet ifakelfake, are the dominant 
producer of events in the diphoton final state and a 
major background for 77 with fake frp- These back- 
grounds come in two different categories; fake pr due 
to energy measurement fluctuations in the calorime- 
ter and fake fr due to pathologies such as picking the 
wrong vertex in events where the true collision did 
not create a vertex or tri-photon events with a lost 
photon. 

To estimate the background due to energy measure- 
ment fluctuations we use the Met Resolution Model. 
The Met Resolution Model considers the clustered and 
unclustered energy in the event and calculates a prob- 
ability, P(|i^'"^* > fJr), for fluctuations in the energy 
measurement to produce equivalent to or larger 

than the measured pr- This probability is then used 

to define MetSig = -log^^o [Pf^fiucy^) ■ Events 

with true and fake ^ of the same value should have, 
on average, different MetSig. For each data event we 
throw 10 pseudo-experiments to generate a ^ and 
calculate its significance, according to the jets and un- 
derlying event configuration. Then we count the num- 



ber of events in the pseudo-experiments that pass our 
MetSig and other kinematic requirements. This num- 
ber, divided by the number of pseudo-experiments, 
gives us the Met Model prediction for a sample. The 
systematic uncertainty on the number of events above 
a MetSig cut is evaluated by comparing the Met Model 
predictions with the default set of model parameters 
to predictions obtained with the parameters deviated 
by ±cr. The total uncertainty is estimated by adding 
the statistical uncertainty on the number of pseudo- 
experiments passing the cuts and these systematic un- 
certainties in quadrature. 

A source of QCD background that is unaccounted 
for by the Met Model is diphoton candidate events 
with event reconstruction pathologies such as a wrong 
choice of the primary interaction vertex or tri-photon 
events with a lost photon. To obtain the prediction for 
the number of events with significant reconstruction 
pathologies in the QCD background at the same time, 
we model the kinematics and event reconstruction us- 
ing a MC simulation of events with in the detector us- 
ing PYTHIA [1^ and a GEANT-based detector simula- 
tion Tg]. We simulate a sample of SM 77 events, with 
large statistics, and normalize to the number of events 
in the presample to take into account jet backgrounds 
which should have similar detector response. Then we 
subtract off the expectations for energy measurement 
fluctuations in the MC to avoid double counting. The 
remaining prediction is due to pathologies alone. The 
systematic uncertainties on this background predic- 
tion include the uncertainty on the scale factor and 
the uncertainty due to MC-data differences in the un- 
clustered energy parameterization and the jet energy 
scale. 

Electroweak processes involving W's and Z's are 
the most common source of real and significant fJr in 
pp collisions. We estimate the background rate from 
decays into both charged and neutral leptons using a 
combination of data and MC methods. There are four 
ways we can get a 77 -I- frp signature in electroweak 
events that decay into one or more leptons: 1) from 
W^77 and 2^77 events where both photons are real; 
2) from Wj and Zj events with a fake photon; 3) from 
W and Z events where both photon candidates are 
fake photons; and 4) tt production and decay. To 
estimate the contribution from the electroweak back- 
grounds we use the Baur [l^l and pythia MC's along 
with a detector simulation, according to their produc- 
tion cross section and k-factors (the ratio of the next- 
to-leading order (NLO) cross section to the leading 
order cross section) , but normalized to data. To min- 
imize the dependence of our predictions on potential 
"MC-data" differences, we normalize, using the rate 
of the number of 67 events observed in the data that 
also pass all signal kinematic cuts, to the number of 
events observed in MC. This ej sample is derived from 
diphoton trigger datasets and the events are required 
to pass the preselection requirements where electrons 
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are required to pass photon-like ID requirements [T3| . 
The uncertainty on the electroweak backgrounds are 
dominated by the 67 normahzation factor uncertainty. 
This includes data and MC statistical uncertainties as 
well as differences in MC modeling. The total uncer- 
tainties also include the MC statistical uncertainties 
and uncertainties on the normalization factors added 
in quadrature. 

Non-collision backgrounds coming from cosmic rays 
and beam-related effects can produce ^^-\-f}r candi- 
dates [§]. These are estimated using the data. Using 
the inclusive 77 sample selection requirements, but re- 
quiring one of the photons to have t^>2b ns we iden- 
tify a cosmic-enhanced sample. Similarly, we utilize a 
beam-related background enhanced sample. We esti- 
mate the number of these events in the signal region 
using the ratio of events outside the timing require- 
ments to events inside the signal region and the mea- 
sured efficiencies of the non-collision rejection require- 
ments [1^. The uncertainties on both non-collision 
background estimates are dominated by the statisti- 
cal uncertainty on the number of identified events. 

After estimating the MetSig distributions for all the 
backgrounds, where the QCD is normalized to the 
data in the low MetSig region where the EWK back- 
grounds are expected to be negligible, the expected 
MetSig distribution for the presample is shown in Fig- 
ured] With these tools in hand we are set to estimate 
the backgrounds for a large variety of kinematic re- 
quirements and move to an estimation of the accep- 
tance for GMSB models in the signal region for use in 
optimization. 



ulation. For the purpose of this analysis we con- 
sider a GMSB model with parameters fixed on the 
minimal-GMSB Snowmass slope constraint (SPS 8) 
that is commonly used 0, @] and take the messenger 
mass scale A/i„=2A, tan(/3)=15, //>0 and the num- 
ber of messenger fields N^—\. The G mass factor 
and the supersymmetry breaking scale A are allowed 
to vary independently. All SUSY production pro- 
cesses are simulated to maximize our sensitivity to 
the model [Hi. 

Since we estimate the sensitivity of the search to 
be equal to the expected 95% C.L. cross section lim- 
its with the no signal hypothesis, we need the un- 
certainties for the luminosity, background and accep- 
tance. The systematic uncertainty on the luminosity 
is taken to be 6% with major contributions from the 
uncertainties on the CLC acceptance from the preci- 
sion of the detector simulation and the event gen- 
erator [1^. The background uncertainty is evaluated 
for every set of cuts in the optimization procedure. 
The systematic uncertainty on the signal acceptance 
for an example GMSB point of m{j(^) = 140 GeV and 
r(x]') ^ 1 ns is estimated to be 6.9% with major con- 
tributions from diphoton ID and isolation efficiency 
(5.4%) and ISR/FSR (3.9%). The uncertainty on the 
NLO production cross section is dominated by the un- 
certainty from parton distribution functions (7.6%) 
and the renormalization scale (2.6%) for a total of 
8.0%. All uncertainties are included in the final cross 
section limit calculation, and we take the acceptance 
and production cross section uncertainties in quadra- 
ture for a total uncertainty of 10.6%. 



YY+Pt analysis in GI\^SB 



CDF Run il Preiiminary, 2.6 fb 

I ' ' ' \ ' ' ' 1^ 



data 

QCD with fake 
EWK with real 
Non-collision 




MetSig 



Figure 1: The background predictions of MetSig for the 
presample. The highest MetSig bin includes all overflow 
events. 



4. GMSB Signal Monte Carlo and 
Systematic Uncertainties 

To estimate the acceptance for GMSB we use the 
PYTHIA event generator as well as a full detector sim- 



5. Optimization and Results 

Now that the background is estimated and the sig- 
nal acceptance is available for a variety of selection re- 
quirements, an optimization procedure can be readily 
employed to find the optimal selection requirements 
before unblinding the signal region. We optimize for 
the following kinematic requirements: MetSig, Ht, 
and A0(7i,72). 

As described in earlier section, the MetSig cut gets 
rid of most of the QCD background with fake pr- 
The Ht cut separates between the high Ex, light final 
state particles produced by GMSB events via cascade 
decays and SM backgrounds, dominated by QCD and 
electroweak backgrounds, which do not have lots of 
high Et objects. The A0(7i, 72) cut gets rid of events 
where two photons are back to back since electroweak 
backgrounds with large Ht are typically a high Et 
photon recoiling against W ev, which means the 
gauge boson decay is highly boosted. Also the high 
Et diphoton with large Ht from QCD background are 
mostly back-to-back with fake fjr or wrong vertex. 

By estimating our sensitivity using the 95% C.L. 
expected cross section limits on GMSB models in the 
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YY+?^ analysis in GMSB 
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Figure 2: The expected 95% C.L. cross section limit as a function of the Ht (a) requirement for a GMSB example point 
(m(xi) = 140 GeV and r(xi) 1 ns). All other cuts held at their optimized values. The optimal cut is where the 
expected cross section is minimized. Indicated in green is the 8.0% uncertainty-band for the production cross section 
and in yellow is the RMS. The N-1 predicted kinematic distribution along with data after the optimized requirements 
are shown in Figure (b). There is no evidence for new physics and the data is well modeled by backgrounds alone. 



no-signal assumption, we find the optimal set of cuts 
before unblinding the signal region. We use the stan- 
dard CDF cross section limit calculator [l^ to cal- 
culate the limits, taking into account the predicted 
number of background events, the acceptance, the lu- 
minosity and their systematic uncertainties. 

For each GMSB point the minimum expected cross 
section limit defines our set of optimal requirements 
for the mass and lifetime combination. The exclusion 
region is defined by the region where the production 
cross section is above the 95% C.L. cross section limit. 
The mass/lifetime limit is where the two cross. Fig- 
ure [21-(a) shows the expected cross section limit as 
a function of a kinematic selection requirement after 
keeping all other requirements fixed at the already 
optimized values, showing it is at the minimum for a 
mass-lifetime combination of m{x\) — 140 GeV and 
t(xi) ^ 1 ns, which is near the exclusion region limit. 

We decided to use a single set of optimal require- 
ments before we open the box based on the obser- 
vation that they will yield the largest expected ex- 
clusion region. We chose: MetSig>3, 77t>200 GeV, 
A0(7i, 72)<7r — 0.35 rad. With these requirements we 
predict a total of 1.38±0.44 background events. The 
dominant electroweak contributions are — > vv^ 
and — !■ ^/i7 which produce a total of 0.26±0.08 
and 0.19±0.10 events respectively. The QCD back- 
ground is dominated by energy measurement fluctu- 
ations in the fr, estimated using the Met Model, to 
have a rate of 0.46±0.24 events. The non-collision 
backgrounds are dominated by cosmic ray which have 
a rate of 0.00llg;[5g^ events. 

After all optimal cuts we open the box and observe 
no events, consistent with the expectation of 1.2±0.4 
events. We show the kinematic distributions for the 
background and signal expectations along with the 



data in Figure [2]- (b). There is no distribution that 
hints at an excess and the data appears to be well 
modeled by the background prediction alone. 

We show the predicted and observed cross section 
limits along with the NLO production cross section, 
which is calculated by multiplying the PYTHIA LO 
cross section calculation by k-factor [23| in Figure [31 
Since the number of observed events is below expecta- 
tions the observed limits are slightly better than the 
expected limits. The Xi mass reach, based on the 
predicted (observed) number of events is 141 GeV/c^ 
(149 GeV/c^), at a Hfetime below 2 ns We show 
the 95% C.L. NLO exclusion region as a function 
of mass and lifetime of Xi using the fixed choice of 
cuts from the optimization for both for the predicted 
and observed number of background events in Fig- 
ure 0]- (a). These limits extend the reach beyond the 
CDF delayed photon results Q and well beyond those 

of D0 searches at r~o <C [tI and the limit from 

Xi 

ALEPH/LEP 4], and are currently the world's best. 



6. Conclusions and Prospects for the 
future 

We have set limits on GMSB models using the 
"f"f + frr final state. Candidate events were selected 
based on 13 times more data, the new frp resolu- 
tion model technique, the EMTiming system and a 
full optimization procedure. We found events using 
2.6 of data in run II which is consistent with 

the background estimate of 1.2±0.4 events from the 
Standard Model expectations. We showed exclusion 
regions and set limits on GMSB models with a Xi 
mass reach of 149 GeV/c^ at a Xi lifetime much less 
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Figure 3: The predicted and observed cross section limits as a function of the Xi mass at a lifetime much less than 1 ns 
(a) and as a function of the Xi lifetime at a mass of 140 GeV/c^ (b). Indicated in green is the 8.0% uncertainty-band 
for the production cross section, in yellow the RMS variation in the expected on the cross section limit. 
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Figure 4: The predicted and observed exclusion region along with the limit from ALEPH/LEP [l] and the 7 + fhr + jet 
delayed photon analysis [![. We have a mass reach of 141 GeV/c^ (predicted) and 149 GeV/c'^ (observed) at the lifetime 
up to 1 ns. The green shaded band shows the parameter space where 0.5 < < 1.5 keV/c^, favored in cosmologically 
consistent models [sj (a). The projected sensitivity to GMSB models with more data. The black dashed line shows the 
prediction of the exclusion region limit after a scaling of the background prediction and the uncertainties for a luminosity 
of 10 fb~^. The blue dashed lines show the prediction of the exclusion region limits from the delayed photon analysis for 
a luminosity of 2 fb~^ and 10 fb~^ respectively taken from Ref. ,8]. (b) 



than 1 ns. Our results extend the world sensitivity to 
these models. 

To investigate the prospects of a search at higher 
luminosity we calculate the cross section limits as- 
suming all backgrounds scale linearly with luminosity 
while their uncertainty fractions remain constant. By 
the end of Run II, with an integrated luminosity of 
10 fb~^, we estimate a mass reach of ~ 160 GeV/c^ 
at a lifetime much less than 1 ns, as shown in Figure ID- 



(b). For higher lifetimes (above ~2 ns) the next gen- 
eration delayed photon analysis will extend the sensi- 
tivity taken from Ref. [1] and then will combine these 
results for completeness. 
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